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Abstract 

We address the study of the classical Gaudin spin model from the bi-Hamil- 
tonian point of view. We describe in details the sl(2) three particle case 

1 Introduction 

The Gaudin models are among the best known integrable lattice models 0. 
Physically, they consist of a system of N s^(2)-valued "spins" A4, possibly coupled 
to an external magnetic field whose strength is assumed to depend on the lattice 
site. Its phase space M can be identified with the n-fold product of sl(2), and the 
Hamiltonian is 



JV N 
H ° = 2 TT ( A i A i) + Tl "( a ^^)' ff > A i G sZ (2) • (LI) 

On M there is a natural Hamiltonian structure, the direct product of the Poisson- 
Lie bracket on sl(2). Complete integrability can be deduced by considering the 
Lax representation of such a model. One considers the Lax matrix L ra t(X) = 

a + £2=1 a^7> and Proves IH that 

a) The Hamiltonian flow associated with Hq is a Lax flow; 

b) integrals of the motion are provided by the spectral invariants of L; 

c) the Lax equations admits an i?-matrix formulation, so that the integrals of the 
motion are in mutual involution. 

In this note we want to address the problem from the standpoint of bi-Ham- 
iltonian geometry, namely to frame these models within the so-called Gel'fand- 
Zakharevich (GZ) |Sj set-up for bi-Hamiltonian integrable systems. Essentially, 
the core of such an approach relies in the extensive use of the Lenard recursion 
relation, associated with a bi-Hamiltonian structure, to generate mutually com- 
muting constants of the motion, via the so-called method of the Casimirs of the 
Poisson pencil. 
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To equip the Gaudin model with a bi-Hamiltonian structure we will use a 
suitable strategy, to be discussed in Section which consists in considering, along 
with the Lax matrix L ra t a polynomial Lax matrix L, po i = rL_i(A — ai)L ra t, and 
recalling that on the space of polynomial matrices a family of mutually compatible 
Poisson structures is defined. The natural Poisson structure we referred above 
turns out to be a specific element of such a family, and so we can define a second 
Poisson structure for the Gaudin model simply choosing another element of such 
a family. Successively, we use this bi-Hamiltonian structure to recover complete 
integrability of the model following the GZ recursion procedure. Finally, we will 
briefly address the problem of Separation of variables, in the bi-Hamiltonian set-up 
discussed injSj. 

Such an analysis can be performed for an arbitrary number of spins (or "par- 
ticles"), and, mutatis mutandis letting the "spins" be sl(n),n > 3 valued. A full 
account of this program is outside the size of this paper. In this note we will use the 
three particle si (2) model for a case study. The above mentioned generalization 
of this approach will be published elsewhere. 



2 The model and the bi-Hamiltonian struc- 
ture 

We consider the three-particle sl(2) Gaudin model, whose degrees of freedom are 
encoded in the Lax matrix with spectral parameter A 



L ra t = 0~ + ^ 



A; 



i=l 



A - en 



where the a, are distinct constant parameters and 



A: 




1,2,3, a 



1 
-1 



(2.1) 



The phase space M is 1 the Cartesian product of three copies of sZ(2); on M one 
can define "standard" Poisson brackets, simply taking the Cartesian product of 
the Lie-Poisson brackets on sl(2): 



hereinafter we will not distinguish between sl(2) and its dual, assuming implicitly to identify the 
two by means of the Killing form. 
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Here and in the sequel we will extensively use a "matrix" representation of these 
brackets (rather, of the Poisson tensor P associated with these brackets), which 
can be explained as follows. 

We identify the tangent and cotangent bundles to M with M itself, so that a 



tangent vector will be written as the vector X 

dF dF OF 



A-t , A 2 , A 3 ) and the differential 



of a function F as dF 



, and notice that the Hamiltonian vec- 

dF, 



dAx dA 2 dA 3 

tor fields associated with the Lie-Poisson brackets are given by A{ = L4j, 
1..3. We can recast in matrix form such a formula as: 



dAi 



X 



A 2 

\A 3 j 



PdF 



I 



\ 



[Ai, 







[A 2 , 







[A3, 



\ 


( dF \ 




dAi 




dF 




dA 2 


/ 


dF 



(2.2) 



Similarly, we will write in an analogous matrix form the linear Poisson brackets 
we are going to discuss in the sequel. 

To endow M with a bi-Hamiltonian structure consider, along with the rational 
Lax matrix L ra t the polynomial Lax matrix 



L 



poly 



i=l 



a,j)L 



rat 



AV + 5 2 A 2 + £iA + £ . 



Explicitly, this change of coordinates on M is given by: 

B = a 2 a 3 Ai + aia 3 A 2 + aia 2 A 3 - s 3 a 

Bi = -(a 2 + a 3 )Ai - (ai + a 3 )A 2 - (a x + a 2 )A 3 + s 2 a 

B 2 = A l +A 2 + A 3 - Sl a 



with 



si 



(ai + a 2 + a 3 ), s 2 = a\a 2 + aia 3 + a 2 a 3 , s 3 = a±a 2 a 3 . 



(2.3) 



(2.4) 



(2.5) 



The rationale for this "coordinate change" is that, on the space of polynomial 
pencils of matrices a family of mutually compatible Poisson brackets are defined 9 
via R-matrix theory. In a nutshell, this family can be described (in our three parti- 
cle case) by saying that there is a map (p from the space of degree three polynomials 
in the variable A to the set of Poisson structures on the manifold of polynomial 
Lax matrices of the form (|2.3j) which sends the monomials A , . . . , A 3 into four 
fundamental Poisson brackets, IL,z = 0, ...,3. These fundamental brackets are 
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written as follows: 





' [Bi,.] 


[S 2 , 


•] 




\ 


n = 


[S 2 ,.] 


[a,. 


] 



















/ 




( ° 




[B 


•] 





n 2 := 


-[B , 


■] - 


[Bi 


•] 







\ 












n. 



n. 



/ 



[S , 










[B2, 



[So, 





k,-] 


[So,-] 
[Si,.] 



/ 



\ [s ,.] [Si,.] [b 2> .] y 



A straightforward computation shows that under the map connecting the Ai's with 
the Bj's, the standard Poisson bracket (|2.2j) is sent into the combination: 



n 3 - S1II2 + s 2 iii - s 3 n 



(2.6) 



of the fundamental brackets. So we can regard P as being associated by <j> with 
the polynomial 

p(X) := (A - ai)(A - a 2 )(A - a 3 ) = A 3 - S1A 2 + s 2 A - s 3 

and we consider the Poisson tensor Q = II2 — S1II1 + S2II0 = <p(q(X)), where 

q(\) = (p(\)/\) + = \ 2 - Sl \ + s 2 . 

Translating back in the coordinates associated with the matrices Aj the tensor Q, 
we get that its components are given by expressions of the form: 



Q 



1 -j 



+ dij [a, 



(2.7) 



\k=l 



where the cj 5 ■ and the dij are somewhat complicated rational functions of the 
parameters 01,02,03, whose explicit expressions are irrelevant here. 

Summing up, we have equipped the phase space M = {sl{2)f of the Gaudin 
model with the Poisson pencil 

P\ = Q — as. 

Quite clearly, the specific choice of the polynomial g(A) is somewhat arbitrary, and 
different choices could be considered. However, in this paper we will stick to this 
choice, and, in the next Section we will study the GZ geometry of such a pencil. 
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3 The GZ analysis of the Poisson pencil 



The Gel'fand-Zakharevich method is a "recipe" to associate with a Poisson pencil 
Q — XP an integrable system (that is, a family of mutually commuting vector 
fields). It is particularly suited for studying Poisson pencils of non-maximal rank. 
Its core can be described as follows: 

Let M be a In + k dimensional bi-Hamiltonian manifold, and let us suppose 
that the rank of P and Q equals 2n. One starts fixing a basis {-f^}i=i...,fc of 
the Casimir functions of P, that is functions satisfying PdHf = 0. Applying the 
second tensor Q to one of such functions one gets (in general) a non trivial vector 
field Xj = QdHf , and tries to find another function Hj- such that QdHf = PdHf . 
If such an equation can be solved for H}, then one generates a new vector field 
applying Q to dHf, and so on and so forth. Supposing one can iterate this process, 
one finds, for each independent Casimir of P, a Lenard chain of vector fields 
X". As a consequence of the bi-Hamiltonian recursion relations, all functions 
Hf commute (w.r.t. both Poisson brackets), even if they do not pertain to the 
same Lenard chain. Complete integrability (of every vector field of the chains) is 
recovered whenever n + k elements of the family of functions H 3 k are functionally 
independent, noticing that k elements of such "fundamental" Hamiltonians are 
provided by the k Casimirs Hf. Indeed, all vector fields of the Lenard chains are 
tangent to the generic symplectic leaf of P, and, when restricted to such a leaf, 
the remaining n functions of the fundamental Hamiltonians provide the required 
n mutually commuting integrals of the motion. 

To apply these ideas to the three particle Gaudin model, we first notice that 
dimM = 9 and the rank both of P and Q equals 6. Actually, a basis Hf of 
Casimirs is given by the three functions 

Hf = TrA 2 t = res A=ai (A - a t )Tr(L 2 rat ) 

Furthermore, applying Q to such functions we get three (independent) vector fields 
X{. Finally, a long but straightforward computation shows that it is possible to 
find three additional functions H\ such that: 

PdHl = Xi, QdH} = 

In other words, we can arrange these six functions in three Lenard chains of the 
form 

Hf H} 

y \ /\ 

Xi 

The new Hamiltonians Hi are a linear combinations of the Casimirs Hf and of 



5 



the other independent spectral invariants Ki of the Lax matrix, given by 
Ki = res A=(l! Tr(0 = 2 ( £ + a^aA,) 



For instance one has 



1 ~~ ? ^27 ^2 K 6 - 1 ) 

ax[ai - a 2 ) z (ai - a 3 ) z 



with 



51 = - 3s 3 ),ti = aia 2 a 3 (ai - a 2 )(ai - a 3 ). 

(i^2 an( i ^3 are obtained by cyclic permutations of the indexes (1,2,3)). Notice 
that there is an invertible linear relation (with constant coefficients) between the 
basis of spectral invariants {Hf,Ki} and the basis of GZ Hamiltonians {Hf,Hf} 
Summing up the bi-Hamiltonian structure we consider gives rise to three Lenard 
chains of "length" one each (that is each one comprising one independent vector 
field). Since the "physical" Gaudin Hamiltonian (Jl.lj) is a linear combination of 
the Ki, H g = ^Ylf=i a i^i an d hence of the GZ Hamiltonians, we have recov- 
ered complete integrability of the model. Finally, notice that we can collect the 
Hamiltonians {H^,H^} into the three polynomials 

F i (\) = \H? + Hl 

and algebraically represent the short Lenard chain(s) depicted above by means of 
the formula (Q - XP) dF^X) = 0. 



4 On Separation of Variables 

In this last Section we will briefly address the problem of Separation of Variables 
for the inhomogeneous Gaudin models. Separation of variables (SoV) for this 
model was proven [T2l PHA |3J |I] , who defined separated variables as coordinates 
of the poles of a suitably normalized Baker-Akhiezer function associated with the 
Lax matrix (|2,1|) . A geometrical approach to the SoV problem, based on bi-Ham- 
iltonian geometry, has quite recently been discussed in the literature jl 11 HI 12*1 15]. 
Here we want to show that SoV for the Gaudin models falls within this scheme. 

The basic geometrical object underlying the bi-Hamiltonian approach to SoV 
is a so-called wJV manifold, that is, the datum of a symplectic manifold (V, lu) 
endowed with an additional (1, l)-tensor N with vanishing torsion satisfying a 
suitable compatibility condition with respect to the symplectic 2-form to (a Nijen- 
huis tensor for short) The content of the "bi-Hamiltonian" SoV theorem (see JH]) 
can be stated as follows: 
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I) On (suitable open sets of) of the ujN manifold V a family distinguished 
canonical coordinates (called Darboux-Nijenhuis (DN) coordinates) is intrinsically 
defined through the spectral properties of its Nijenhuis tensor N jlUj . 

II) A Liouville integrable system, characterized by a complete set of mutually 
(w.r.t. the Poisson tensor w -1 ) integrals of the motion {-Hj}i=i ... iv is separable in 
DN coordinates if and only if these Hamiltonians commute also w.r.t. the second 
bracket defined by 

Q = N-uj- 1 . (4.1) 

As we have seen, the Gaudin models (as it is the case for a number of integrable sys- 
tems related with R-matrix theory and/or reductions of soliton equations) admit 
a (reasonably) natural bi-Hamiltonian formulation on a bi-Hamiltonian manifold 
where none of the brackets of the pencil is non degenerate. So, the first problem 
to be tackled in such an instance is the so-called reduction problem, that is to 
concoct from the (degenerate) bi-Hamiltonian structure of the problem a structure 
of uN manifold on suitable submanifolds of M, in such a way to satisfy condition 
II above. 

To show that this can be successfully achieved in the Gaudin models we will 
follow a recipe discussed in [21^], specifically suited for GZ systems. We consider 
the Poisson pencil Q — XP, and fix a basis {H®, . . . , H?} of Casimirs of P, where 
we assume that k = corank(P) = corank{Q). The symplectic leaves S c of P 
are thus characterized by the equations Hf = const j, and come equipped with a 
natural symplectic structure. It holds [3] the 

Proposition: In the above geometrical set-up, let Z\ , . . . , Zh be a family of vector 
fields on M, transversal to the symplectic leaves of P such that the functions van- 
ishing along the distribution < Z{ > generated by the Zi 's are a Poisson subalgebra 
both w.r.t. P and Q. Then the generic symplectic leaf of P has the structure of an 
u)N manifold. The Nijenhuis tensor N is defined, via Eq. 1)4 by the restriction 
of the following (modified) Poisson tensor Q = Q — Yli=i Q(d{Hj,° P )) A Zi, where 
Z{ are a normalized basis (possibly defined on an open subset of S c ) for < Zi >, 
i.e. Zj(Hf) = 5ij. Condition II above is automatically satisfied. □ 

We will now apply this scheme to the three particle sl(2) Gaudin model, en- 
dowed with the bi-Hamiltonian structure Q — XP introduced in Section [2j As we 
have seen, the GZ structure of the problem is quite simple. We have three Casimir 
polynomials Fi(X) of degree one, (meaning that, indeed corank(P) = corank(Q) = 
3) and so we have to find a three-dimensional distribution satisfying the properties 
of the Proposition recalled above. 

The idea to solve this problem is very simple, and relies on the following obser- 
vation on the (ordinary) Lie- Poisson brackets on a single copy of si (2). With the 

notations of Section |2] the Poisson bracket of two functions F, G on M, is given 

dF dG OF dG 

by -f-F, G] = Tr(— — • L4, — — -1) = —Tr(A ■ \-r-r, -r—r})- Parametrizing the generic 
y dA 1 dA u K l dA dA u 
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element A of sl(2) as A 



h/2 f \ Q 
we consider the vector field Z = — . 

ae 



e -/i/2 

We notice that differentials of functions vanishing along Z admit a very simple 

( 

matrix representation. Indeed Z is matricially represented as Z(A) = 

V 1 

fdF\ OF 
and so Z(F) = iff I I = 0, i.e., iff — — - lies in the lower Borel subalgebra 

V / 1 2 

b of sl{2). 

Let i 7 , G be functions such that Z{F) = Z{G) = 0, and let us compute 
Z({F,G}). Thanks to the Leibniz property of the Lie derivative and the fact 
that Z is a constant vector field we have that 

f) F f)G 

Z({F,G})=-Tr(Z(A) ■[—,—]) (4.2) 

which vanishes as well since b_ is indeed a Lie subalgebra of sl{2). Finally we 
notice that the Casimir C of the Lie-Poisson bracket is given by C = h? /2 + 2ef, 
so Z{C) = 2f and hence Z is generically transversal to the symplectic foliation 
of the Lie Poisson brackets; in particular, the normalized generator Z is given by 
1/2/Z. 

To use this simple result in the case of the three particle sl(2) Gaudin model we 

d 

consider the vector fields Zi = — — ,i = 1, 2, 3.. The differentials of functions F van- 

dF OF dF 

ishing along < Zi > are represented by triple of matrices dF = (77^—, -7—^ , 77^—) 

oA\ oA 2 0A3 

dF 

with 7—— G b_, i = 1, 2, 3. The fact that such functions are a Poisson subalgebra 
dAi 

for P is self evident. To ascertain that the same is true for Q one simply has notice 
that, using the explicit expressions l|2.7j) . the brackets {F, G}q =< dF, QdG > are 
given by the multiple sum 

Since a is a constant and the vector fields Zi admit the matrix representation 
Zi(Aj) = 5ijZ(A) we see that the Lie derivatives Zi({F, G}q) is a multiple sum 
of terms like those of Eq. (|4.2|) . and so vanish whenever Z(F) = Z(G) = 0. This 
proves that on the symplectic leaves of P the bi-Hamiltonian pencil Q — XP induces 
an uN structure whose DN coordinates separate the Hamilton-Jacobi equations 
of the Gaudin models. 

Acknowledgments: This work is an outgrowth of a long standing research 
project of one of us (GF) with F. Magri and M. Pedroni. 
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